Electronic and structural characterization of divacancies in irradiated graphene 
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We provide a thorough study of a carbon divacancy, a fundamental but almost unexplored point 
defect in graphene. Low temperature scanning tunneling microscopy (STM) imaging of irradiated 
graphene on different substrates enabled us to identify a common two- fold symmetry point defect. 
Our first principles calculations reveal that the structure of this type of defect accommodates two 
adjacent missing atoms in a rearranged atomic network formed by two pentagons and one octagon, 
with no dangling bonds. Scanning tunneling spectroscopy (STS) measurements on divacancies 
generated in nearly ideal graphene show an electronic spectrum dominated by an empty-states 
resonance, which is ascribed to a spin-degenerated nearly flat band of 7r-electron nature. While the 
calculated electronic structure rules out the formation of a magnetic moment around the divacancy, 
the generation of an electronic resonance near the Fermi level, reveals divacancies as key point 
defects for tuning electron transport properties in graphene systems. 

PACS numbers: 68.37.Ef, 73.20.Hb, 73.22.Pr, 75.70.Rf 



Graphene is a unique material from an electronic point 
of view. The ultrarelativistic nature of its charge carri- 
ers [l[ and the robustness of its electronic coherence Q 
make it the ideal candidate in the forthcoming era for na- 
noelectronics. However, all these remarkable electronic 
and transport properties are subjected to the presence of 
ubiquitous disorder. By instance, extrinsic/intrinsic de- 
fects commonly present in graphene are considered the 
limiting factor for electronic transport through charged 
impurities 0, Q , rippling [Ej or resonant scatterers • 
These defects lead to substantial changes in the topol- 
ogy of its low-energy electronic bands, which are at the 
origin of many of its unique properties. Taking advan- 
tage of the key role of defects in low dimensional carbon 
systems, a new route based on defect-engineering is be- 
ing developed to broaden the functionalities of graphene 
d, 0| . Adatoms, vacancies and Stone- Wales (SW) de- 
fects are the most common defects invoked to act as 
building blocks in structurally tailored graphene. Both 
vacancy-type and SW defects can be artificially created 
in graphene by electron or ion irradiation and visualized 
with atomic resolution by high-resolution transmission 
electron microscopy [To| and scanning probe microscopy 
[ill [l2| . In particular, STM has demonstrated to be a 
powerful tool to explore the local electronic structure of 
such point defects at the atomic scale [13l-ll5]. 



The potential effects of topological defects in carbon- 
based systems have been comprehensively studied by the- 
ory in the last years, predicting already new functional 
properties like the opening of an electronic gap in the 
Dirac bands in graphene with SW defects |18l or aiming 
to explain the unexpected magnetic ordering [16|, LL7j ob- 



served in defective carbon systems [19l-l25j. Therefore, a 



detailed picture of all types of point defects existing in 



graphene is mandatory in order to understand their par- 
ticular impact on the materials properties and to be able 
to selectively tailor them. From the experimental point 
of view, although STM is an ideal technique to study the 
properties of point defects in graphene systems as single 
entities, few works are yet reported, mainly restricted to 
single atomic vacancies. 

In this Letter, we present a combined experimental and 
theoretical characterization of a divacancy, a point defect 
expected to be stable and common in ion irradiated car- 
bon samples [26|, |27j . Divacancies are also expected to be 
of fundamental relevance regarding the electron trans- 
port properties of these systems [28| being considered, 
for instance, to be the key defect in order to tune the 
conductance of carbon nanotubes [29|.The experimen- 
tal data here reported were acquired at 6K by using a 
home-made low temperature STM (LT-STM) in ultra- 
high- vacuum (UHV) conditions [30|. The local density 
of states (LDOS) of the sample was studied with atomic 
precision by measuring differential conductance (dl/dV) 
spectra using the lock-in technique (f = 2.3 kHz, Vmod = 
1.5 mV). All STM/STS data were acquired and processed 
using the WSxM software [3l[ . In this work we have per- 
formed various density functional [32| calculations (DFT) 
using the SIESTA code (33l, l34j which uses localized or- 
bitals as basis functions |35j. We use a double ( polar- 
ized basis set, non-local norm conserving pseudopoten- 
tials and a functional for the exchange and correlation 
including van der Waals interactions [H, [I?} • The calcu- 
lations are performed with stringent criteria in the elec- 
tronic structure convergence (down to 10 -5 in the density 
matrix), Brillouin zone sampling (up to 900 k-points), 
real space grid (energy cut-off of 500 Ryd) and equilib- 
rium geometry (residual forces lower than 0.02 eV/A). 
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FIG. 1. (color online) Divacancy in four different graphene 
systems: (a) lML-G/SiC(0001) (V=+20 mV), (b) 2ML- 
G/SiC(0001), (V=+500 mV), (c) HOPG (V=+500 mV) and 
(d) 4-5ML-G/SiC(0001) ( V=+190 mV) . (a)-(d) Size = 6x6 
nm 2 . All bias voltages refer to the sample, (e) STM image 
(26 x 26 nm 2 ) showing two divacancies in 4-5ML-G/SiC(000T) 
(V=+280 mV). The direction of graphene lattice vectors ai 
of the outermost layer is indicated at the top right part of the 
image. 



We deliberately generated divacancies in several 
graphene systems by means of Ar + irradiation in UHV 
(140eV), following the same experimental procedure as 
in references [HI, Hif . We chose four different graphene 
systems: monolayer and bilayer graphene on SiC (0001), 
4-5 layers of graphene grown on SiC(OOOl) and the sur- 
face of UHV-exfoliated highly oriented pyrolitic graphite 
(HOPG). The methods used to grow epitaxial graphene 
on both SiC orientations (Si and C face) are described in 
[38l . [39! ] and [!o| respectively. After the irradiation pro- 
cedure we imaged the samples by LT-STM to explore the 
created defects. Among the artificially generated defects, 
we were able to identify a common point structure exist- 
ing in all these graphene surfaces, which is displayed in 
Figs. l(a)-(d). According to the STM images, the most 
characteristic feature of such defect is a two-fold sym- 
metry involving a complex electronic pattern. Another 
common characteristic of all these defects is the orienta- 
tion of the mirror plane parallel to their long axis which is 
rotated 30 with respect to the graphene lattice vectors 



(see fig 1(e)). The similar shape of the defects and their 
equal orientation with respect to the honeycomb lattice, 
strongly suggest that they all belong to the same kind 
of vacancy- type defect. Interestingly, the occurrence of 
this type of defect in each one of the graphene systems 
studied here turns out to be quite different. While for the 
SiC systems this kind of defect was frequently observed 
after the irradiation procedure, in the case of HOPG sur- 
faces such defect was extremely rare, indeed, the defect 
shown in Fig. 1(c) was the only one found on the HOPG 
surface after extensive measurements in many different 
irradiated HOPG surfaces. The reasons explaining this 
strong unbalance in the divacancy population remain un- 
clear to us and they are most likely related to the different 
underlying environment for each graphene system. 

In order to further investigate this kind of defect we 
focused on graphene grown on the C terminated face of 
SiC, where the rotational disorder of the graphene lay- 
ers electronically decouple bands for a large interval of 
rotational angles [4ll443|. This can lead to a stacking of 
undistorted and nearly isolated graphene sheets in this 
system, as shown by the experimental observation of ideal 
non-doped graphene cones around K points [3]. There- 
fore, surface STM measurements on multilayer graphene 
on SiC (0001) can be considered as a model experimental 
approximation to the study of ideal neutral monolayer 
graphene. Figure 1(e) shows a large scale STM image 
of an irradiated surface of 4-5ML-graphene/SiC(000T) 
with two of such vacancies present in the terrace (la- 
beled as A and B). In addition to the atomic graphene 
lattice, two other super-periodicities can be observed in 
the STM topograph. Such modulations, known as moire 
patterns [45| , are frequently found in graphene grown on 
SiC (0001) surfaces [40| and they are due to the rotational 
misalignment of the upmost graphene layers. The moires 
observed here arise from two different rotations of 20 
and 5 between graphene layers from an AA stacking 
leading to a periodicity of 7.1Aand 29. 5 A, respectively. 
The larger moire is compatible with a commensurate su- 
perstructure (n,m) = (6,7) with a moire vector V= nai 
ma2, ai being the graphene lattice vectors (see inset fig 
le) and n,m integers. 

The electronic decoupling due to the rotation between 
graphene layers is illustrated in figure 2. First, we have 
calculated a bilayer graphene structure forming a (2,3) 
commensurate moire using the SIESTA code, as indi- 
cated above. Figure 2(a) shows a simulated STM im- 
age, within the Tersoff-Hamann approximation (4^, of 
the relaxed structure, in which the moire can be readily 
observed. The calculated moire corrugation (5.3 pm) re- 
sults from both topographical (3.2 pm) and electronic 
(2.1 pm) contributions which are in phase. The cor- 
responding low-energy electronic band structure of this 
moire bilayer is represented in fig. 2(b), together with 
the calculated ones for the monolayer and AB bilayer. 
While the linear band dispersion of monolayer graphene 
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FIG. 2. (color online) (a) Simulated STM image of a commensurate moire (2,3) formed by two graphene layers rotated G 
= 13.17° from an AA stacking. The unit cell (yellow) has a periodicity of 10. 7A. Parameters: 5x5nm 2 , V = +200mV. (b) 
Electronic band structure calculated for monolayer graphene (dash blue line), bilayer AB (green line) and the (2,3) moire 
(red dot line), (c) Corresponding DOS for the calculated structures, d) Typical dl/dV spectrum acquired on the surface of 
4-5ML-Graphene/SiC(0001) at 6K. Plots showed in (c) and (d) have the same energy range than in figure (b). 



becomes quadratic in the vicinity of the Dirac point (Ed) 
in a Bernal AB bilayer, for such energies the electronic 
structure of graphene is restored in the bilayer when both 
layers are rotated, as already reported (4ll-l43|. Interest- 
ingly, although the calculated DOS shows in all cases a 
similar "V" shape in this energy range (see fig. 2(c)), only 
the rotated bilayer shares a vanishing DOS at Ed with 
the monolayer graphene. Our tunneling spectroscopy 
experiments performed on different moires in graphene 
on SiC(0001) agree with these calculations. Figure 2(d) 
shows a dl/dV curve measured on a pristine region. The 
typical V-shaped and the vanishing differential conduc- 
tance at Ed (in this case at -20mV with respect to 
the Fermi energy (E F ) due to residual n-doping of our 
SiC(0001) samples) confirms an effective decoupling of 
the graphene surface layer with respect to the substrate, 
allowing a direct comparison of the experimental data 
with calculated structures for a single monolayer. 

The experimental results presented in figure 1 do not 
reveal unambiguously the atomic structure of the stud- 
ied defect. The main contribution to the STM images 
around defects is electronic so they only provide indirect 
fingerprints of the atomic rearrangement. However, the 
structure of the defect likely involves vacancy sites as ion 
irradiation on graphene systems mainly produces voids in 
the atomic network. More precisely, recent ion bombard- 
ment simulations on graphene performed by Lehtinen et 
al. , predict that low-energetic impacts should mostly pro- 
duce single and double vacancies rather than more com- 
plex structures like Frenkel pairs or SW defects. [27]. The 
calculated ratio of mono- and divacancies in graphene, 
according to our irradiation parameters (E(Ar + ) = 140 
eV at 20° off normal), is roughly 2:1 [47J. As stated pre- 
viously, single atomic vacancies are the only well studied 
point defects in graphene systems by STM, showing a 
characteristic electronic pattern with a well defined three- 
fold symmetry (lll-[l3|. in contrast to the present defect. 
Therefore, we performed DFT calculations focused on 
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FIG. 3. (color online) (a) Relaxed structure of the calculated 
divacancy where all sigma bonds saturate leading to the (585) 
structure. Dashed lines show the mirror planes in the defect, 
(b) 2.57 x 2.57 nm 2 unit cell used for the calculation. The 
right panel shows the comparison between experimental [(c) 
and (d)] and calculated [(e) and (f)] STM images at both 
sample bias polarities. All images have a size of 3.3 x 3.3 



divacancies structures to complete the characterization 
of the two-fold symmetry vacancy-type defect here ob- 
served. There are various possible atomic configurations 
saturating all dangling bonds from a honeycomb lattice 
where two carbon atoms are removed. The resulting lat- 
tices are locally formed by non hexagonal polygons that 
saturate all a dangling bonds. The simplest of these 
atomic structures for a reconstructed divacancy is the 
(585) defect, shown in figure 3(a). The rearranged pla- 
nar structure encloses a central octagon and two opposing 
fivefold rings. We have chosen this configuration as the 
initial structure for the calculation due to its two-fold 
symmetry and the relative orientation of the two orthog- 
onal mirror planes to the honeycomb lattice, ° and 30 °, 
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FIG. 4. (color online) Consecutive tunneling differential con- 
ductance spectra acquired on the divacancies labeled as A 
and B in fig. 1(e). A reference dl/dV curve taken on the bare 
surface is also shown. 



as measured in the STM images. Indeed, previous DFT 
calculations of this divacancy in the LDA approximation 
disp lay similarities with the present experimental images 
[22j| . In our monolayer supercell DFT calculation, and 
in order to minimize the interaction between defects in 
neighboring cells, we have considered large skewed unit 
cells. In particular, we present the results for the unit cell 
defined by the basis vectors bi= 5ai and b 2 = 7a 2 (see 
Fig. 3(b)). Similar results have been obtained with other 
choices of the unit cell provided they are large enough. 
As can be noticed, the simulated images obtained from 
the relaxed structure (Figs. 3(e)-(f)) capture the main 
features of the experimental STM images (Figs. 3(c)-(d)) 
at both polarities: a central bright lobe surrounded by a 
complex electronic (v // 3xa/3)R30 ° pattern with twofold 
symmetry in an equal orientation. This excellent agree- 
ment allows us to unambiguously identify this defect with 
a (585) divacancy. 

Topological defects in graphene systems break the 
translational symmetry of the crystal leading to profound 
changes in the low-energy electronic bands. Thus, to 
complete the characterization of this class of divacancies, 
we have investigated their local electronic structure by 
means of STS experiments and DFT calculations. Re- 
cent theoretical works 2o|, HH EE nave predicted the 
existence of quasilocalized states in the vicinity of the 
Ep for different types of topological defects in graphene 
such atomic vacancies, point impurities or edges. Indeed, 
these electronic modifications are expected to be at the 
heart of the striking behavior observed in functionalized 
graphene regarding its transport and magnetic proper- 
ties. These states have been experimentally observed 
by some of us in single atomic vacancies in HOPG and 
graphene on Pt(lll) surfaces [HI EH- Our STS findings 
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FIG. 5. (color online) (a) Calculated band structures and 
the corresponding DOS of pristine monolayer graphene (dots) 
and monolayer graphene with a (585) divacancy (solid), (b) 
Comparison between the experimental STS curves shown in 
fig. 4 for the divacancy B with the calculated DOS around 
the E F . 



in the present case are summarized by the local tunnel- 
ing spectra taken at the center of the A and B divacan- 
cies on the SiC(0001) region shown in fig. 1(e). Figure 
4 shows dl/dV curves acquired consecutively on both A 
(continuous red) and B (point blue) defects within 0.5 eV 
around Ep together with a reference curve measured with 
the same tip on a neighboring pristine graphene region. 
While the latter spectrum shows a "V" shape with a DOS 
vanishing at Ep, w Ep characteristic of ideal graphene as 
discussed above, it is clear that the presence of vacancies 
profoundly alter the LDOS. dl/dV spectra recorded on 
A and B divacancies show a resonance located in empty 
states and centered at +150mV, which present a rich in- 
ternal structure. Although both defects share the main 
features in their electronic spectra, and thus the LDOS 
reflects qualitatively an equal electronic structure, some 
subtle differences in the internal structure of the reso- 
nances were systematically distinguished. The calculated 
band structure for a divacancy in a monolayer with a 
(5,7) unit cell is shown in Fig. 5(a) along with the band 
structure of the defect free monolayer. In this figure the 
effect of the presence of the divacancy is apparent; it 
induces an almost dispersionless 7r-character band inter- 
acting with the continuum, giving rise to a 0.3 eV width 
asymmetric resonance in the DOS (see Fig. 5(a)). Thus, 
as can be noted in Figure 5(b), the calculated total DOS 
for the (585) divacancy is in nice agreement with the ex- 
perimental dl/dV curves. However, although the width 
of the resonances are similar in both cases, the internal 
structure of the experimental resonance is not clearly re- 
flected in the calculated DOS. Highly demanding further 
calculations including a second rotated graphene layer 
and involving a very large number of carbon atoms could 
be required in order to understand this subtle effect. Fi- 
nally, although we have performed a spin-resolved calcu- 
lation, the obtained charge distribution is identical for 
both spin states, leading therefore to a non magnetic so- 
lution. 

In summary, we present a comprehensive characteri- 
zation of the geometrical and electronic structure of a 
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(585) divacancy, a common point defect in irradiated 
graphene systems. This defect has a planar and non- 
hexagonal rearranged structure with no dangling bonds 
and thus it is expected to present a low chemical reactiv- 
ity compared to other vacancy- type defects. Diferential 
conductance spectra taken on divacancies generated at 
the surface of 4-5ML- graphene/SiC(0001) reveal the ex- 
istence of an electronic resonance centered at +150mV. 
The calculated band structure for the divacancy in a sin- 
gle layer of graphene fully supports the STS observations, 
proving the existence of a slightly dispersive electronic 
state of 7r nature at this energy. Indeed, according to 
our DFT calculations, this type of divacancy does not 
show a magnetic character. The resonance associated to 
(585) divacancy at low energies is expected to limit the 
electron mobility [28], being thus an excellent candidate 
to be used to functionalize graphene by ion irradiation 
methods. 
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